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ABSTRACT 

The  effect  of  surface  tension  on  the  shape  of  a  two-dimensional  jet 
emerging  from  an  orifice  is  considered.  It  is  shown  that  the  slope  of  the 
surface  profile  of  the  jet  is  not  continuous  at  the  separation  points.  Both 
velocity  and  curvature  are  infinite  at  these  points.  The  problem  is  solved 
numerically  by  series  truncation.  Jet  profiles  are  presented  for  various 
values  of  the  surface  tension.  In  addition,  perturbation  solutions  for  both 
small  and  large  values  of  the  surface  tension  are  derived,  jr-— 
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SIGNIFICANCE  AND  EXPLANATION 

The  classical  Kirchhof f 1  solution  for  the  shape  of  a  two-dimensional  jet 

«eerging  from  an  orifice  yields  infinite  curvature  at  the  separation  points. 

It  is  shown  that  this  singularity  is  not  removed  by  including  surface  tension 

in  the  boundary  condition.  On  the  contrary  surface  tension  makes  the  problem 

more  singular  by  introducing  a  discontinuity  in  the  slope.  This  result  is  in 

4  5 

agreement  with  Vanden-Broeck's  '  findings. 

The  problem  is  solved  numerically  by  collocation  for  arbitrary  values  of 
the  surface  tension.  In  addition  perturbation  solutions  for  both  small  and 
large  values  of  the  surface  tension  are  presented. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  NEC,  and  not  with  the  author  of  this  report. 


THE  EFFECT  OF  SURFACE  TENSION  ON  THE  SHAPE  OF  THE  KIRCH HOFF  JET 


Jean-Marc  Vanden-Broeck* 

1.  INTRODUCTION. 

The  classical  Kirchhof f 1  solution  for  the  shape  of  a  two-dimensional  jet 
emerging  from  an  orifice  yields  infinite  curvature  at  the  separation  points. 

A  similar  singularity  arises  in  cavitating  flow  problems  such  as  those 
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considered  by  Ackerberg  ,  Cumberbatch  and  Norbury  and  Vanden-Broeck  '  . 

4 

Vanden-Broeck  considered  the  influence  of  surface  tension  on  the 

cavitating  flow  past  a  flat  plate.  He  provided  numerical  and  analytical 

evidence  that  the  slope  of  the  free  surface  is  discontinuous  at  the  separation 

points.  Thus  the  inclusion  of  surface  tension  did  not  remove  the  infinite 

curvature  at  the  separation  points.  On  the  contrary  it  made  the  problem  more 

singular  by  introducing  a  discontinuity  in  slope.  These  results  were 

generalized  by  Vanden-Broeck^  for  the  cavitating  flow  past  a  curved  obstacle. 

In  the  present  paper  we  investigate  the  effect  of  surface  tension  on  the 

Kirchhof f1  jet*  A  numerical  scheme  based  on  series  truncation  is  presented  to 

solve  the  problem  for  arbitrary  values  of  the  surface  tension.  The  numerical 
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results  are  qualitatively  similar  to  those  obtained  by  Vanden-Broeck  '  .  The 
slope  is  not  continuous  at  the  separation  points.  Both  velocity  and  curvature 
are  Infinite  at  these  points.  In  addition  perturbation  solutions  for  both 
small  and  large  values  of  the  surface  tension  are  presented. 
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The  probla  is  formulated  In  Motion  2.  The  numerical  results  ore 
presented  in  Section  3.  The  perturbation  calculations  ere  derived  in  Section 


4  end  5. 


2.  rowroiATiow. 


We  consider  e  two-dimensional  jet  asarging  from  an  orifice  (see  Figure 
D*  The  orifice  is  assumed  to  be  hole  in  a  plane  wall  of  small  thickness. 

The  effects  of  gravity,  viscosity  and  compressibility  are  neglected.  Far 
downstream  from  the  orifice  the  speed  in  the  interior  of  the  jet  is  equal  to  a 
constant  D  and  the  two  streamlines  bounding  the  jet  are  straight  and 
parallel.  We  denote  by  L  the  thickness  of  the  jet  far  downstream  from  the 
orifice. 

We  define  dimensionless  variables  by  choosing  L  as  the  unit  length 
and  U  as  the  unit  velocity.  We  introduce  the  potential  function  4  and  the 
stream  function  4.  Without  loss  of  generality  we  choose  4*0  at  the 
separation  points  and  4*0  on  the  streamline  IJ.  It  follows  from  the 
choice  of  the  dimensionless  variables  that  t  ■  1  on  the  streamline  IJ*. 

We  denote  the  complex  velocity  by  u  -  iv  and  we  define  the  function 
T  -  i8  by  the  relation 


T-iO 


u  -  iv  -  a  (1) 

We  shall  seek  x  -  i8  as  an  analytic  function  of  f  ■  4  +  in  the  half 
plane  4  <  0.  The  complex  potential  plane  is  sketched  in  Figure  2. 

On  the  surface  of  the  jet  the  Bernoulli  equation  and  the  pressure  jump 
due  to  surface  tension  yield 


12  T  _  1  2 

i"  *;* ■s° 


(2) 


Here  q  is  the  flow  speed,  K  the  curvature  of  the  jet  surface,  T  the 
surface  tension  and  P  the  density.  Zn  dimensionless  variables  this  becomes 


Hare  a  is  the  Weber  number  defined  by 

OLU2 

ft  m 

T 


The  kinesuitic  condition 

in 

IS  and 

IS*  yield 

0  - 

0, 

♦  -  0 

♦  <  o 

(6) 

0  “ 

0, 

♦  -  1 

♦  <  o 

(7) 

This  completes  the  formulation  of  determining  the  function  t  -  10.  For 
each  a,  x  -  10  must  be  analytic  is  the  half  plane  ♦  <  0  and  satisfy  the 
boundary  conditions  (3 ) ,  (4),  (5)  and  (6). 

3.  NUMERICAL  RESULTS 

Following  Birkhoff  and  Zarantonello6  we  define  the  new  variable  t  by 
the  relation 

l(t.l)  ,a, 

This  transformation  naps  the  flow  domain  onto  the  unit  circle  in  the  complex 
t-plane  so  that  the  walls  go  onto  the  real  diameter,  and  the  free  surfaces 
onto  the  circumference  (see  Figure  3). 

we  introduce  the  function  Q(t)  by  the  relation 

x  -  10  -  in  t  -  jp-  Jlnd  -  t2>  +  Q(t>  (9) 

Here  0  is  the  value  of  0  at  ♦  ■  ♦  •  0.  The  function  ft(t)  is  bounded 
and  continuous  in  the  unit  circle  It |  <  1,  and  analytic  in  the  interior. 

The  conditions  (6)  and  (7)  show  that  Q(t)  can  be  expressed  in  the  form  of  a 


Taylor  expansion  in  even  powers  of  t.  Hence 


T  -  10  *  tn  t  -  ~  in(1  -  t  )  +  l  U  .t 
w  *•_  n+1 

n*0 


The  function  (10)  satisfy  (5)  and  (6).  The  coefficients  Un  and  the  constant 

0  have  to  be  determined  to  satisfy  (3).  The  condition  (4)  will  then  be 

automatically  satisfied  by  symmetry. 

■  i  i  o 

We  use  the  notation  t  *■  |t|e  so  that  points  on  SJ  are  given  by 
t  *  ei0,  0  <  a  <  Using  (8)  we  rewrite  (3)  in  the  form 

_  cos  0  T  d8  a  .  2t  .... 

HiTo  *  1  <*  *  ”  <11> 

Here  t(o)  and  9(0)  denote  the  values  of  r  and  0  on  the  free  surface 


We  solve  the  problem  approximately  by  truncating  the  infinite  series  in 

(10)  after  N  terms.  We  find  the  N  coefficients  U  and  the  constant  0 

n 

by  collocation.  Substituting  t  -  e  °  into  (10)  and  taking  the  real  and 
imaginary  parts  we  obtain 


-  ~  log (2  sin  o)  +  ][  U  ..cos  2no 


-  ?  0  *  3  +  l  D„+18in  2n° 

n»0 


We  now  introduce  the  N+1  mesh  points 

*  2(d  +  i)  ^  ”  2^'  ^  ^  ^ 

A* 

**  <30 

Using  (12)  and  (13)  we  obtain  [T(o)]  and  [7s’]  in  terms  of  the 

ao  <x*0 j 

coefficients  UR  and  the  constant  0.  Substituting  these  expressions  into 
(11)  at  the  points  oT  we  obtain  N+1  nonlinear  algebraic  equations  for 


the  N  +  1  unknowns  0,  ,  n  *  1,...,N  and  0.  We  solve  this  system  by 
Newton's  method. 

Once  this  system  is  solved  for  a  given  value  of  a,  the  shape  of  the  jet 
is  obtained  by  numerically  integrating  the  exact  relations 

”7  ■  7  tgoe  Tcos  0  (15) 

<30  n 

a/ 

~  tgoe  Tsin  0  (16) 

Here  x(o)  and  y(o)  are  the  values  of  x  and  y  on  the  free  surface 
SJ.  Relations  (15)  and  (16)  follow  directly  from  (1)  and  (8). 

Some  of  the  numerical  results  are  shown  in  Table  1.  These  values  are 
correct  to  the  number  of  decimal  places  shown.  As  n  increases,  the 
coefficients  0R  decrease  rapidly.  Typical  profiles  are  shown  in  Figure  4. 

Our  numerical  results  (see  Table  1)  indicate  that  0*0  for  all  values 
of  a  *  Therefore  the  slope  of  the  surface  profile  of  the  jet  is 
discontinuous  at  the  separation  points  for  all  values  of  a  *  -.  Both 
velocity  and  curvature  are  infinite  at  these  points. 

In  Figure  5  we  present  a  graph  of  0  versus  a.  As  a  varies  between 
0  and  -,  0  varies  continuously  from  ~  to  0. 

For  a  -  -,  the  condition  (11)  reduces  to  the  free-streamline  condition 
x  *  0.  The  solution  is  then 

0-0,  Un  -  0,  n  -  1,2,3,...  .  (17) 

Substituting  (17)  into  (10)  we  obtain  the  following  exact  solution 

rm  -  i©^  -  in  t,  0-0  (18) 

This  is  the  classical  Kirchhoff^'s  solution.  Here  the  subscript  •  denotes 
a  - 

As  a  ♦  0,  the  profile  of  the  jet  approaches  two  horisontal  straight 
lines.  The  solution  is  then 


ft  -  J  '  O,  -  in  2 


U  *  0  f  n  B  2^ 3| •  •  •  • 
n 


Substituting  (19)  into  (10)  we  obtain  the  following  exact  solution 


*.  -  ■  *”  ^2  <»•> 

Here  the  subscript  0  denotes  a  *  0. 

The  contraction  ratio  C  of  the  jet  is  defined  as  the  ratio  of  the  width 
of  the  jet  far  away  from  the  orifice  to  the  width  of  the  orifice. 

For  a  ■  •  we  obtain  from  (18) 

c-  “  7TT  ~  °*61 

For  a  •  0  we  obtain  from  (20) 

c0  "  1  * 

In  Figure  6  we  present  numerical  values  of  C  versus  a.  As  a  decreases 
from  infinity,  the  contraction  ratio  increases  sonotonically  from  to  CQ. 


4.  PFHTURBATIOH  SOUTTIOH  FOR  tt  LARGE 

For  o  -  •,  the  solution  is  given  by  (1.8).  Using  (8)  and  (18)  we 


obtain  after 


algebra 


♦)”1/2  as  ♦  ♦  0,  t  -  0  (21) 

Therefore  the  Kirchhoff1  solution  yields  infinite  curvature  at  the  separation 
points . 

A  perturbation  solution  for  a  large  can  be  derived  by  using  the  method 

of  matched  asymptotic  expansions.  The  details  of  the  calculation  follow 

4  5 

closely  the  work  of  Vanden-Broeck  '  .  Therefore  they  will  not  be  repeated 
here. 


-12- 


The  asymptotic  results  for  a  largo  confirm  that  tho  slop#  is  not 
continuous  at  tho  ooparation  points.  In  particular  tho  asynptotic  expression 
for  0  is 

-1/2 

0  ~  *(2o)  as  a  ♦  •  (22) 

Halation  (22)  is  shown  in  Figure  5.  The  numerical  results  are  in  good 
agreement  with  (22)  for  a  large.  For  a  -  20  the  value  of  0  predicted  by 
(22)  agrees  with  the  numerical  results  within  two  percent. 


5.  PKRTPRBATIOW  SOLPTIOH  FOR  a  SMALL 

For  o  »  0,  the  solution  is  given  by  (20).  He  seek  a  solution  for  a 
small  in  the  form  of  an  expansion  in  powers  of  a.  Thus  we  write 


Here  end 

we  have 


T  *  T0  +  aTi 

e  -  e0  ♦  oe, 

are  defined  by  (20) 


do  2w 


♦  0(e2)  (23) 

♦  0(o2)  (24) 

.  Substituting  (23)  and  (24)  into  (11) 

cos  a  (25) 


integrating  (25)  and  using  the  condition  we  obtain 


w 

2 


~  sin  o-l  .  w 

0  - - - - a  +  - 

1  2r  2 


(26) 


Relation  (26)  and  the  definition  of  0  imply 


0  -  -  -  2- 
p  2  2w 


(27) 


Relation  (27)  is  shown  in  Figure  5.  For  a  -  1  the  value  of  0  predicted  by 
(15)  agrees  with  the  nusterical  results  within  two  percent. 
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